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Abstract. The Antennae is a pair of late type spirals in 
the course of merging. The interaction triggered an on- 
going strong burst of star formation that also produced 
a large number of Young Star Clusters (YSCs), many of 
which seem to be young Globular Clusters (GCs). The ob- 
served Luminosity Function of the YSC system is a power- 
law. We use evolutionary synthesis models for star clusters 
in comparison with HST WFPC observations of the YSC 
system in the Antennae to analyse the mass function of 
this ~ 2 ■ 10 8 yr young star cluster system. Properly ac- 
counting for age spread effects by individually age-dating 
the YSCs we find that the intrinsic Mass Function (MF) 
is log-normal in shape and in intriguing agreement with 
the MF of old GC systems. We discuss this MF in the 
context of cluster formation and dynamical effects in the 
tidal field of the parent galaxy and speculate about its 
future evolution. 



1. Introduction 

The Antennae (NGC 4038/39 = Arp 244) is a pair of late 
type spirals in the course of merging. It is the youngest in 
Toomre's (1977) dynamical age sequence of 11 interacting 
and merging systems and comfortably nearby. This makes 
it a well-studied system all over the wavelength range from 
X-ray through radio wavelengths. The interaction has trig- 
gered a strong burst of star formation that has also pro- 
duced a large number of ~ 700 bright star clusters as 
first observed with HST by Whitmore & Schweizer (1995) 
(WS95). Many of those bright clusters seem to be young 
Globular Clusters (GCs) due to their small effective radii 
and high luminosities. The Luminosity Function (LF) of 
this bright star cluster system looks like a power law 
$(L) ~ L-i-78±o.05 with no hint of a turnover at fainter 
magnitudes down to the completeness limit of V = 22.3 
mag which, at the distance of The Antennae (19.2 Mpc 
for H = 75), corresponds to M v = -9-6 mag (WS95). In 
Fritze - v. Alvensleben (1998) (hereafter Pap. I) I anal- 
ysed WS95's HST data with evolutionary synthesis mod- 
els for single burst single metallicity populations (SSPs) 
using stellar evolutionary tracks for various metallicities 



from the Geneva group and a Scalo IMF from 0.15 to 
60 M Q (Fritze - v. Alvensleben & Burkert 1995, FB95). 
In analogy to the YSC system in NGC 7252 which, like 
the Antennae, is an Sc - Sc merger, we assume that the 
metallicity of the YSC population that forms out of the 
spirals' ISM is ~ i • Z Q (Fritze - v. Alvensleben & Ger- 
hard 1994). In Pap. I individual ages are determined from 
(V — I) colours for the 550 star clusters with I - band 
detections. The age distribution clearly reveals two popu- 
lations of clusters, ~ 70 old GCs from the original spirals, 
and ~ 480 Young Star Clusters (YSCs) with ages in the 
range (0 — 4) • 10 8 yr formed in the ongoing interaction- 
induced starburst. Only the secondary population of YSCs 
will be considered in the following. Meurer (1995) was the 
first to point out the possible importance of age spread 
effects on the future evolution of the LF of a YSC sys- 
tem. With individual YSC ages and the fading in vari- 
ous passbands as given by our SSP models we were able 
to calculate the future evolution over a Hubble time of 
the LF and of the colour distribution of the YSC sys- 
tem under the unrealistic assumption that all clusters will 
survive. At an age of ~ 12 Gyr, when age differences of 
the order of 10 8 yr do no longer play any role, the LF 
is shown to be undistinguishable from a typical GC LF 
(i.e. Gaussian with a turnover at My ~ —6.9 mag - ap- 
propriate for the metallicity [Fc/H] ~ —0.3 (Ashman et 
al. 1995) - and cr(Mv) = 1.3 mag). The colour distribu- 
tion will by then also be compatible with the one observed 
on old GC systems. While accounting for stellar mass loss, 
this modelling, however, did not take into account dynami- 
cal effects like two-body relaxation, dynamical friction, or 
disk shocking that might act heavily on a YSC system 
over a Hubble time. Already for the Milky Way, where 
the potential is rather well known, it is difficult and not 
uncontroversial to model the dynamical processes on indi- 
vidual clusters (e.g. Chernoff & Weinberg 1990, Fukushige 
k Heggie 1995, Gnedin & Ostriker 1997, Vesperini 1997). 
This is, of course, even more difficult, if not impossible, in 
an ongoing merger like The Antennae. Before one can try 
and quantify the efficiency of various destruction mecha- 
nisms, the intrinsic MF underlying the presently observed 
LF has to be determined. Knowing individual cluster ages 
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offers the possibility to use M/L values from evolutionary 
synthesis models to derive the Mass Function (MF) un- 
derlying the presently observed LF and this is what we 
attempt in this Letter. Since the age distribution of the 
YSCs is strongly peaked within < 2 • 10 8 yr and the YSCs 
have a mean distance to the galaxy center of <~ 3.5 kpc, 
a YSC on average cannot have had more than 1 or 2 rev- 
olutions. We do not expect the MF therefore to already 
be significantly affected by cluster destruction processes. 
Rather we expect the presently derived MF to reflect the 
MF produced by the cluster formation process. 

The mass spectra of molecular clouds, molecular cloud 
cores, open clusters, and the LF of giant HII regions (in 
non-interacting galaxies) all are power laws with expo- 
nents a in the range a ~ —1.5 ... —1.7 (e.g. Solomon et al. 
1987, Lada et al. 1991, Kennicutt 1989, see Harris & Pu- 
dritz 1994 and Elmcgrccn & Efrcmov 1997 for overviews) 
as is the mass spectrum of open clusters in the Milky Way 
and the LMC (e.g. van den Bcrgh & Lafontaine 1984, El- 
son & Fall 1985). Both the LF and the MF of old GC sys- 
tems are Gaussians with typical parameters (My) ~ —7.3 
mag, (7 <~ 1.3 mag, and (Log (M/M Q )} ~ 5.5, a <~ 0.5, re- 
spectively (e.g. Ashman et al. 1995). 

The question immediately arises: Is the transition from 
a power law molecular cloud mass spectrum to a Gaus- 
sian old GC mass spectrum performed in the star/cluster 
formation process or by secular destruction effects within 
a GC system ? Or, else, is already the mass spectrum 
of molecular clouds or molecular cloud cores different in 
strongly interacting and starbursting galaxies from what 
it is in normal spirals ? 



2. The Mass Function of the YSCs in The 
Antennae 

On the basis of individual YSC ages we use our SSP mod- 
els giving M/L a in the passbands A ~ UBVRIK as a func- 
tion of time to derive masses of individual clusters from 
their observed V - luminosities. 

This is done for all the 393 YSCs with ages < 4 • 10 8 
yr and V - luminosities brighter than the completeness 
limit My = —9.6 mag. It is stressed that our model M/L 
- values include the decrease in cluster mass due to stellar 
mass loss (cf. Pap. I), but not that due to the evapora- 
tion of stars from the cluster. The MF we recover in this 
way from the presently observed LF is presented in Fig.l. 
A Gaussian with (Log (M Y sc/M©)} = 5.6 and a = 0.46, 
normalised to the number of YSCs in the histogram, is 
overplotted. The intrinsic MF we obtain for the YSCs 
brighter than the completeness limit in The Antennae 
clearly looks log-normal in shape with the maximum at a 
mean YSC mass of ~ 4 • 10 5 M Q . Stellar mass loss within 
the clusters from the present mean age of ~ 2 • 10 8 yr 
through an age of ~ 12 Gyr will lead to a decrease in 
mass of < 15 % for a Scalo - IMF (< 10% for Salpeter). 
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Fig. 1. Mass Distribution of the YSCs in The Antennae: 
393 YSCs brighter than the completeness limit. A Gaus- 
sian with (Log (Mysc/M q )) = 5.6 and a = 0.46 is over- 
plotted, normalised to the number of YSCs in the his- 
togram. 



Thus, without any destruction or evaporation effects 
the mean mass of the secondary GC system in The Anten- 
nae at the age of a Hubble time would be ~ 3.4 • 10 5 M . 
A cluster with this mean mass would have My = —6.9 
mag at an age of 12 Gyr according to our models. This is 
the position of the maximum of the Gaussian YSC LF we 
obtain at a hypothetic YSC age of 12 Gyr (cf. Fig. 6a in 
Pap. I). The agreement is no surprise since this is, in fact, 
the way how we obtained the parameters for the Gaussian 
in Fig.l. We stress that these parameters are derived from 
the LF in Pap. I and are not the result of any fit to the 
cluster MF in Fig.l. 

Remarkably enough, the parameters of this Gaussian 
- which in Fig.l is seen to reasonably describe the MF 
of the secondary GC population - are quite similar to 
those given by Ashman et al. (1995) for the Milky Way 
and M31 GC systems. Using evolutionary synthesis re- 
sults from Worthey (1994) for M/Ly, Ashman et al. find 
(Log (M/Mq)) = 5.47 and a = 0.50 for the Milky Way 
and (Log (M/M Q )) = 5.53 and a = 0.43 for M31. The MF 
in Fig. 1 thus seems compatible with the bulk of the YSCs 
really being young GCs rather than open clusters or asso- 
ciations, as was already indicated by their small effective 
radii and high luminosities (WS95). 

From the model side, uncertainties in the determina- 
tion of YSC ages (and hence masses) on the basis of their 
(V — I) colours are dominated by uncertainties in the YSC 
metallicities. Age uncertainties due to metallicity uncer- 
tainty (i • Z Q < Zysc 5; Z q ) are estimated to be of the 
order of ±15 %. The uncertainty in M/L at Z = i • Z Q 
due to the age uncertainty is ~ 8 % and the uncertainty 
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in M/L at all ages < 4 • 10 8 yr due to the metallicity un- 
certainty is < 5 %. This leads to an overall uncertainty in 
the M/L of - 10 %. 

Measurement uncertainties are ~ ±0.2 mag for the ob- 
served (V — I) colours and < 0.15 mag for V magnitudes 
for YSCs brighter than the completeness limit. The uncer- 
tainties that the inhomogeneous reddening over the body 
of NGC 4038/39 brings along for the derived cluster ages 
and masses are difficult to quantify. Only a global average 
value is given for the internal reddening of the YSC pop- 
ulation in NGC 4038/39 (E V -i ~ 0.3 mag (WS95)) and 
applied before our age-dating, but dust lanes and unre- 
laxed structures are seen all over NGC 4038/39. The very 
good agreement (< 10 7 yr) between ages determined from 
(U — V) and (V — I) colours, however, indicates that for 
the bright clusters seen on the short U exposure the red- 
dening does not seem to significantly deviate from the 
average reddening we use. 

It should be noted in this context that inclusion of 
YSCs fainter than the completeness limit - which tend to 
have significantly larger observational uncertainties - does 
not affect the agreement with the Gaussian in Fig.l or its 
parameters (beyond the normalisation to the number of 
clusters in the histogram). 

We conclude that the secondary GC population 
in The Antennae is formed with a log-normal mass 
distribution very similar to the one in the Milky 
Way or M31. It is not necessarily the secular evolu- 
tion but rather the cluster formation process that 
produces the Gaussian mass spectrum observed on 
old GC systems. 

Since uncertainties in the MF from observational er- 
rors cannot be calculated in any straightforward way in 
the analysis presented here, we are currently trying an 
independent approach. We draw YSCs at random from 
different intrinsic MFs (Gaussians, power-laws), randomly 
assign ages to them from different age distributions (clus- 
ters formed uniformly over the burst duration or at an 
increasing or decreasing rate), and calculate their present 
luminosities and colours to which, then, observational er- 
rors can be added, again at random from the observed 
luminosity-dependent error distribution. Comparison of 
the resulting model LFs and colour distributions with the 
observed ones should then allow to constrain the intrinsic 
MF (Kurth et al. , in prep.). A somewhat similar proce- 
dure is used for YSCs in NGC 7252 and NGC 3921 by 
Miller & Fall (1997), who find that power-laws are pre- 
ferred over Gaussians for the MFs of these YSC systems. 

3. Comparison with YSCs in NGC 7252 and 
NGC 3921 

NGC 7252 and NGC 3921 are the two oldest merger rem- 
nants from Toomre's (1977) list. Large enough YSC sys- 
tems have been detected in both of them to define the 
bright end of their LFs (Whitmore et al. 1993, Miller et 



al. 1997, Schweizer et al. 1996). Distances, however, are 
larger than to NGC 4038/39 by factors 3 and 4 for NGC 
7252 and NGC 3921, respectively, pushing the complete- 
ness limit to significantly higher luminosities. The higher 
mean ages of the YSCs in these cases (650 — 750 Myr for 
NGC 7252 and 250 - 750 Myr for NGC 3921, depending 
on metallicity) add another argument in favour of them 
being young GCs, since they have survived 3> 10 • t cross 
(cf. Schweizer et al. 1996). 

The LF we calculate for the YSCs in The Antennae 
at an age of ~ 1 Gyr does show a marginal turnover at 
the expected (My) — —9.9 mag for YSCs brighter than 
the completeness limit, indicating that by this age the 
distortion of the LF with respect to the MF due to age 
spread effects from the finite YSC formation timespan (of 
the order 200 — 400 Myr) is already less important. 

In order to estimate if a turn-over in their YSC LFs 
could be expected in NGC 7252 and NGC 3921, we assume 
that their YSC systems have the same MF as the YSC sys- 
tem in The Antennae. Then, we can calculate the mean ab- 
solute magnitude (My) from (Mysc) at the above-quoted 
mean ages. We obtain (My) = — 10 . . . — 9.5 mag and 
(M v ) = -9.5 ... - 9 mag for YSCs in NGC 7252 and 
in NGC 3921, respectively. These luminosities are close to 
the 90 % completeness limiting magnitudes of —9.5 (PC) 
and -8.5 (WF) for NGC 7252 (Miller et al. 1997) and of 
-9.0 for NGC 3921 (Schweizer et al. 1996). Together with 
the difficulty of disentangling the old GC population from 
the YSCs at these advanced ages the fact that no turnover 
is detected in the LFs for NGC 7252 and NGC 3921 does 
not seem to rule out a Gaussian MF similar to the one we 
obtain in The Antennae. 

4. Dynamical Evolution of the YSC System 

Quantitatively, nothing is known about the external dy- 
namical effects on GC systems in interacting and merg- 
ing galaxies. For the Milky Way, where dynamical mod- 
elling on the basis of the Galactic potential is possible and 
has been done by many groups (e.g. Chernoff & Weinberg 
1990, Gnedin & Ostriker 1997, Fukushige & Heggie 1995, 
Vesperini & Heggie 1997, Vesperini 1997, 1998) it is clear 
that the GC system observed today is only "the hardiest 
survivors of a larger original population" (Harris 1991). 
Vesperini & Heggie (1997) present N-body simulations in- 
cluding effects of stellar evolution, two-body relaxation, 
disk shocking, and the tidal field of the Milky Way. Study- 
ing the secular evolution of a number of GC systems with 
different initial MFs, Vesperini (1998) shows that if the 
initial MF is log-normal then this log-normal shape 
and its parameters are conserved over a Hubble 
time despite the destruction of a large part (~ 50 %) of 
the original GC population. While evaporation and dy- 
namical friction preferentially destroy low and high mass 
clusters, respectively, both processes balance each other in 
the case of a log-normal (= equilibrium) initial GC MF, so 
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that no net selective destruction of specific GC masses re- 
sults. If the GC destruction processes were similar in the 
Antennae and in the Milky Way, and if, as indicated in 
Fig.l, the initial GC MF really were close to Vesperini's 
equilibrium MF, then the Gaussian MF could survive a 
Hubble time with its parameters (Log (M/M Q )) and a 
essentially unchanged despite the likely destruction of a 
large fraction of the YSC system seen today. 

It will be interesting to analyse more YSC systems in 
the way we did in order to see if and in how far the YSC 
MF is universal or depends on parameters of the progen- 
itor galaxies or the interaction event. Moreover, studying 
secondary GC systems on an age sequence of interacting, 
merging and merger remnant galaxies should allow to di- 
rectly "observe" both the time evolution of the LF and 
of the MF and thus offer a unique possibility to study 
the effects of dynamical processes in situ. The turnover of 
an old GC population at My ^ — 7.2 mag occurs around 
V ~ 24.5 at Virgo cluster distances and is well within the 
reach of 10m telescopes. HST imaging allows to identify 
clusters in regions not too crowded for ground-based MOS. 

5. Discussion and Open Questions 

The uncertainties in the metallicities and - more impor- 
tant in the individual reddening of YSCs in The Anten- 
nae can be substantially reduced with spectroscopic obser- 
vations of YSC selected from the HST images. With MOS 
facilities on 10m class telescopes this should be possible 
in the nearest future. 

One question that poses itself immediately in the con- 
text of the cluster formation mechanism is whether the 
molecular cloud mass spectrum in a strongly interacting 
and starbursting galaxy like The Antennae does really 
have the same power-law form as that in an undisturbed 
spiral forming stars at a level orders of magnitude lower? 
Elmcgrccn & Efremov (1997) pointed out that the high 
pressure of the ambient ISM produced in a strong interac- 
tion might favor the production of more massive clouds. 
Information about the molecular cloud mass spectrum in 
merger-induced strong starbursting systems seems a pre- 
requisite for the study of the star and cluster formation 
processes. 

First spectroscopic observations of a handful of YSCs 
in NGC 7252 (Schwcizcr & Seitzer 1993, 1998) confirms 
the metallicity of (| — 1) • Z (cf. FB95) we predicted on 
the basis of the progenitor galaxies' ISM abundances. This 
enhanced metallicity (with respect to the primary GC 
population) raises the question in how far the secondary 
GC formation process is comparable to the primary one 
in the early Universe? 

6. Preliminary Conclusions 

1. The MF of the YSCs in The Antennae seems to be log- 
normal with parameters (Log (M/M Q )) and a very similar 



to those of the Milky Way GC system. 

2. This suggests that the cluster formation process and 
not the dynamical evolution produce the Gaussian MF. 

3. The close agreement of the YSC MF we obtain with 
Vesperini's equilibrium initial GC MF seems to indicate 
that the shape and parameters of this MF may survive 
a Hubble time despite destruction of a large fraction of 
today's YSCs. 

4. In the older merger remnants NGC 7252 and NGC 3921 
the completeness limit is close to the turn-over luminosity 
expected in case their MFs were similar to the one in The 
Antennae. 

As long as the impact of observational colour and lu- 
minosity uncertainties on the MF we derive cannot be 
quantified rigorously, our conclusions have to remain pre- 
liminary. 
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